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Abstract. The primordial helium abundance Y p is commonly inferred from abundance determinations in low- 
metallicity extragalactic Hn-regions. Such determinations may be subject to systematic uncertainties that are 
investigated here. Particular attention is paid to two effects: ic/-corrections for "imperfect" ionization structure 
leading to significant amounts of (unobservable) neutral helium or hydrogen and "tcf -corrections due to non- 
uniform temperature. Model H n-regions with a large number of parameters are constructed and it is shown 
that required corrections are almost exclusively functions of two physical parameters: the number of helium- 
to hydrogen- ionizing photons in the illuminating continuum Q(He°)/Q(H°), and the ratio of width to radius 
5rs/rs of the Stromgren sphere. For clouds of sufficient helium-ionizing photons Q(He°)/Q(H°) > 0.15 and non- 
negligible width of the Stromgren sphere a significant overestimate of helium abundances may result. Such clouds 
show radiation softness parameters in the range —0.4 < log r\ < 0.3 coincident with the range of r\ in observed 
Hn-regions. Existing data of Hn- regions indeed seem to display a correlation which is consistent with a typical 
~ 2 — 4% overestimate of helium abundances due to these effects. In case such an interpretation prevails, and 
in the absence of other compensating effects, a significant downward revision of Y v may result. It is argued that 
caution should be exercised regarding the validity of commonly quoted error bars on Y p . 
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<*. 1. Introduction 
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The observational determination of the primordial light el- 
ement abundances provides an important key test for the 
validity of the standard model of Big Bang nucleosynthe- 
sis and the cosmic baryon density. Within the last years, 
a large number of high-quality observations of hydrogen- 
and helium- emission lines originating from extragalactic 
giant Hn-regions and compact blue galaxies has led to a 
number of determinations of the primordial helium mass 



It seems presently feasible to determine the primor- 
dial D/H ratio to unprecedented accuracy by observa- 



fraction Y v (e.g. 


3 eimbert & Torres-Pcimbcrt 1974 


; 


<unth 


fe Sargent 1983; 


Pagel et al. 1992; 


[zotov et al. 1994 


1997; 


Olive ct al. 1997 


; Izotov & Thuan 1998; Pcimbcrt et al. 



200C ) . Though quoted statistical errors of these Y p deter- 
minations are often < 1% due to the size of the employed 
samples, it is by now widely accepted that errors in the 
inferred Y p may be dominated by systematic uncertain- 
ties. Only in this way, may one understand how different 
groups arrive at either h igh Y v « 0.244 (Izotov & Thuan 



199S) or low Y p « 0.234 (Olive et al. 1997; Peimbert et al 



200C) values which deviate by several times the quoted 
statistical error. 
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tions of quasar absorption line systems (Buries et al. 1999; 
|Q'Meara et al. 2000| ; |Tytler et al. 2000| ). These observa- 
tions favor a "low" primordial deuterium abundance, i.e. 
D/H w 3 x 10~ 5 , implying a baryonic contribution to the 
critical density of ^h 2 w 0.02 (with h the Hubble con- 
stant in units of 100 km s _1 Mpc _1 ), and a primordial 
helium abundance of Y p « 0.247, within the context of a 
standard Big Bang nucleosynthesis scenario. The accurate 
observational determination of Y p would not only give an 
independent "measurement" of Q,bh 2 , but would also allow 
for a test of the validity of a standard Big Bang nucleosyn- 
thesis scenario. Such an undertaking, nevertheless, would 
require Y p determinations with statistical and systematic 
errors < 1% , a magnitude which requires investigation 
of a large number of possible errors which may enter the 
analysis. 

Beyond the pure random errors that generally underlie 
any measurement, there are several sources of systematic 
uncertainty in commonly used observational techniques 
and subsequent analysis of spectra of giant H n-regions for 
the determination of helium (and metal) abundances. To 
a first approximation, helium abundances in nebulae can 
be easily inferred by observing the relative line fluxes of 
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lines produced during helium recombinations (Pcimbert & 
Torres- Peimbert 1974) and lines produced during hydro- 



gen recombinations (e.g. H/3), 
/(He I, A) Jn e n Hc+ afJ\,T) 



He 



oc 



(1) 



I{U(3) Jn e n H+ a c H s (R(3,T) " H 

when cloud temperature T and the effective recombina- 
tion coefficients af s (cf. Oster brock 1989) are known 



The influence of systematic errors starts at the correction 
for effects of reddening, underlying stellar absorption (i.e. 
correcting for absorption troughs in the stellar continuum 
at the position of the emission lines), and fluorescent en- 
hancement of helium lines (e.g., through absorption of He I 
A3889 by the metastable 2 3 S level of Hei and re-emission 
as He I A7065). The use of Eq. ([!]) also presumes that 
there is no contribution to the observed line radiation from 
collisional excitation, an effect which becomes important 



partici larly at higher densities. Izotov fc Thuan (1998 ) 
have devised a method to correct these potential errors 
simultaneously by employing several He I emission lines in 
the analysis. This method (and others) has also been re- 
cently ^crhkiahyas^ (see 



also Sassclov fc Goldwirth 1995| ), with the result that 
residual errors may still be appreciable. Uncertainties in 
the theoretically computed effective recombination coeffi- 
cients could also introduce systematic errors as large as 
~ 1.5% ( |Benjamin et al. 199S| ). 

The largest systematic errors, nevertheless, may re- 
sult from the idealizing assumptions that Hn-regions are 
clouds at constant temperature and density, with sim- 
ple geometry, and with helium- and hydrogen- Stromgren 
sphere radii coinciding to within one percent. These as- 
sumptions enter implicitly in an analysis which follows 
the strategy of Eq. ([j]) . In this paper the validity of two of 
these assumptions are tested in detail: the assumption of 
constant temperature (the deviation of this case is com- 
monly referred to as the existence of "temperature fluctu- 
ations" ) and the assumption of equal Stromgren spheres 
(correction for this effect is commonly achieved by fac- 
toring in "ionization correction factors (icf)"). The lat- 
ter (icf) effect is usually only considered by estimating 
the quality of the incident continuum by the "Radiation 



Softness Parameter" rj (Vilchez & Pagel 1988; Skillman 
198 ~ 



j,nd based on model calculations by Stasihska (1990) 
it has been assumed that for sufficiently hard radiation 



(log?/ < 0.9) «c/-corrections are negligible (< 1%, Pagel 
et al. ( 1992| ) ). Nevertheless, the potential importance of 



ionized helium in regions where hydrogen is almost com- 
pletely neutral has been recently stressed by a number of 



groups (Armour et al. 1999; Viegas et al. 2000; Ballantyne 
et al. foOOp. (The p roble m was already noted earlier e.g 



by Stasihska (1980| ), Dinerstcin & Shields (1986), and 
Pcna (1986).) New model calculations therein confirmed 
the presence of this problem. The present work is distin- 
guished from prior analyses by the consideration of a much 
larger number of model H n-regions and incident stellar 
spectra. It also reaches a somewhat different conclusion 
about the importance of this problem. 



Effective recombination coefficients a cS employed to 
infer ionic abundances from observed emission line ra- 
tios are functions of the electron temperature within the 
emitting volume. Temperatures in H n-regions are usu- 
ally inferred from flux ratios of collisionally excited oxy- 
gen lines (i.e. [O m]AA4959, 500 7/4363) and are mostly ap- 
proximated to be constant (cf. Peimbert et al. (200C) for 
an analysis going beyond this). Since collisionally excited 
lines are exponentially sensitive to temperature, a temper- 
ature determination by such methods is systematically bi- 
ased to the hottest parts of the H n-region. However, those 
regions which are somewhat colder may contribute most 
of the helium- and hydrogen recombination lines. A sys- 
tematically overestimated electron temperature will also 
lead to an overestimate of Y p . The problem of the choice 
of an appropriate mean electron temperature f or recom- 
bination coefficients was discussed in detail by Peimbert 



(1967 ) . A discussion of the temperature effect pertaining 



to the determination of Y p has also appeared in the work of 
|Steigman et al. (1997 ) , who gauge the effect by somewhat 
arbitrarily changing inferred temperatures. The present 
work analyses this problem by constructing models for 
Hn-regions with a photoionization code (Cloudy 90.05 
( |Fcrland 1997| )). 

The general strategy of this paper is as follows. The 
simplest possible, spherically symmetric, Hn-regions at 
constant (and low) density are constructed and their emis- 
sion spectra are calculated with the help of a photoioniza- 
tion code. These models include a variety of parameters. 
Degeneracies of the resultant emission spectra to these pa- 
rameters are outlined. Ranges of typical emission line ra- 
tios are inferred from observed Hn- regions, and a model 
H II- region is accepted only if it falls into these ranges. 
Correction factors icf and tcf are defined which quantify 
the appropriate ionization correction and correction for 
temperature inhomogeneity needed to infer the "true" he- 
lium abundance. These factors equal unity for the ideal 
ionized, "one zone" cloud. Finally, the obtained icf and 
tcf correction factors will be compared to observational 
diagnostic tools (i.e. Radiation Softness Parameter and 
[O m1A500 7/[Oi1A6300 ratios) as well as to the H II- region 
sample by Izotov k, Thuan (1998| ). 



In Sect. |2| the model calculations and the parameters 
employed in the calculations are described. Detailed re- 
sults for the required ionization correction and tempera- 
ture correction factors are presented in Sect. || and dis- 
cussed more broadly in Sect. ||. Conclusions are drawn in 
Sect. |. 

2. Models of H n-Regions 

The model H n-regions are calculated using the photoion- 
ization code CLOUDY Version 90.05 (Ferland 1997). In all 
models a spherical symmetric distribution of gas around a 
compact, point-like source of ionizing radiation at the cen- 
ter is assumed. All models have a constant (total) hydro- 
gen number density of jih = 10 cm -3 . This relatively small 
density is chosen to ensure that collisional enhancement 
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of helium recombination lines is negligible. Further, it is 
assumed that effects of radiative transfer can be treated in 
the On-The-Spot approximation. Thus the emission of re- 
combination lines may be described using the well-known 
Case B limit. Note that the photoionization code Cloudy 
has been slightly modified to employ helium emissivities 
relative to H/3 as given in Pagel et al. (1992). In the limit 
of small density these are practically identical to those 
given in Benjamin et al. (199E). The model clouds inves- 



tigated here are ionization bounded, i.e. the edge of the 
cloud is denned by the Stromgren sphere where the de- 
gree of ionization of hydrogen drops rapidly to zero. In 
this case all UV-photons with hv > 13.6 eV are absorbed 
and eventually converted into lower energy photons. 

The a dopted abundances of heavy elements are those 
given by Brcsolin et al. (1999 ). These are based on the 



compilation by Grevesse fc Anders (1989 ), but with deple- 
tion on dust grains taken into account. In order to simulate 
H n-regions with metallicities far below solar, all elements 
heavier than helium are scaled down in abundance by the 
same factor. The employed metallicities in the model cal- 
culations range between Z/Z & — 1/36 and 1/12, corre- 
sponding to 10 6 O/H = 23.64 to 70.93. The effects of dust 
grains, themselves, on, for example, the reddening of lines 
or the heating/cooling balance, is not specifically consid- 
ered. Assuming that the density of dust scales with metal- 
licity, these effects should be fairly small for the clouds 
considered. The helium abundance was kept constant for 
all models at a value of He/H = 0.0776. 

Given the hydrogen number density n H , the spectrum 
of the ionizing background, and the metallicity of the 
cloud, there are three remaining parameters describing the 
simulated H n-regions. These are Q(H) the number of pho- 
tons with hv > 13.6 eV radiated by the source in unit time, 
ro the inner radius of the cloud, and e the volume filling 
factor of gas at density n H taking into account the effect 



of gas condensa tions in H n-regions (see e.g. psterbrock 



& Flat her 1959 ). Nevertheless, it will be shown that the 
variation of only one of these parameters (Q(H), ro, or e) 
is necessary to generate Hn- regions with different emis- 
sion spectra. The quantities Q(H), n H , e, and ro define the 
ionization parameter Us, i.e. the approximate number of 
ionizing photons per hydrogen at the hydrogen Stromgren 
sphere rg: 



Us 



9(H) 
47r r|n H c 



(2) 



Here the Stromgren sphere radius r§ may be somewhat ar- 
bitrarily denned as the location where n H o /n H + = 1 . The 
parameter Us fully determines the structure of a spher- 
ically symmetric cloud for a given chemical composition 
and incident continuum. Thus clouds with the same Us, 
but different rt H , Q(H), and e are self-similar in their ion- 
ization structures. Except for the total luminosity, such 
self-similar clouds will result in the emission of essentially 
identical nebular spectra. This assumes that the effects of 
collisional excitation on the heating/cooling balance and 
the line emission are negligible. This is typically applicable 



for densities n H <SC 100cm -3 . Moreover, as long as models 
have r§ > ro, the exact choice of the parameter ro is of 
secondary importance for the calculation of nebular emis- 
sion spectra. This is primarily because the regions in the 
very interior of the cloud contribute very little to the to- 
tal line emission due to the small amount of gas involved. 
When rs 3> ro one can show from the balance of hydrogen 
ionizations to recombinations that Us is proportional to 

U S oc (Qe 2 n H ) 1/3 • (3) 

Note that this definition of the ionization parameter leads 
to an increase of Us with increasing rt H in contrast to 
the ionization parameter at the inner radius of the cloud 
Uq = rt 7 /rt H (ro), which is often used as the input pa- 
rameter to specify the relative photon to gas density in 
plane-parallel clouds. Figure [l] illustrates this proportion- 
ality Eq. (|^) of Us for a sample of models with varying 
Q(H), r , and e. The few models that do not satisfy Eq. (|3|> 
are those where rs ~ ro, in particular, where shell geome- 
try pertains. In what follows, it is therefore not necessary 
to vary Q(H), ro, and e independently, but rather analyze 
clouds with different Us, the physical parameter deter- 
mining the emission spectra when the geometry of the 
H n-regions is characterized by a sphere. 



Fig. 1. Relationship between 
Q, e, and Us for a sam- 
pl e of models ionized by 
a Mihalas (1972) continuum 
with T cB = 45 000 K. The 
few models that are more 



17 



O 

3> 16 



15 



-4 -3 
log U s 



distant from the line of pro- 
portionality are those where 
the assumption of a solid 
sphere with ro <C rs is no 
longer valid (see text). 



Incident Continuum: The effects studied here are partic- 
ularly sensitive to the shape of the incident continuum 
at energies higher than 1 Ryd. Spectra with a larger con- 
tribution of helium ionizing photons cause an ionization 
structure where the Stromgren sphere of helium is equal or 
even larger than the Stromgren sphere of hydrogen. These 
are spectra of O and early B stars with effective temper- 
atures above T eff ~ 40 000 K. The effective temperature of 
an entire star cluster has these high values during the early 
stages of its evolution. In this case the spectrum is domi- 
nated by the most massive and, thus, hottest stars. Since 
massive stars burn out more quickly than less massive 
ones, at later times, the main contribution to the ionizing 
spectrum is increasingly provided by cooler stars. This, in 
turn, implies a radiation spectrum of the cluster which is 
softer. In this paper several different spectra for the ioniz- 
ing radiation are used. These include two different spectra 
for single stars at various effective temperatures, in partic- 
ular, spectra computed from the non-LTE plane-parallel 
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stellar atmospheres by Mihalas (1972) and the LTE plane- 
parallel atmospheric grids by Kurucz (1991). These con- 
tinua are shown in Figure |^ for one effective temperature. 
For reference, a blackbody spectrum of the same temper- 
ature is also shown in the diagram. All plotted spectra are 
normalized to have the same number of ionizing photons. 
From left to right, the vertical lines indicate the threshold 
energies for ionization of H°, He , and He + , respectively. 
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Fig. 2. Single star spectra by Mihalas (1972) and Kurucz 
(199l|)] Both are shown for an effective stellar tempera- 
ture T eft = 45 000 K and for the same number of hydro- 



gen ionizing photons Q(H) = 10 



51 , 



For comparison, a 



blackbody spectrum with the same parameters is shown 
as well. Fluxes are given at a distance of 73 pc to the star. 
From left to right, the vertical lines show the threshold 
energies for ionization of H°, He and He + . 
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Fig. 3. The continua emitted by stellar clusters of dif- 
ferent ages with metallicity Z/Z & = 1/20 (see text for 
details) . The dotted line shows a blackbody emission with 
the same effective temperature while the vertical dotted 
lines indicate the ionization frequencies of H , He and 
He+. 



ionizing photons (hv > 1.8 Ryd) to that of hydrogen ion- 
izing photons (hv > lRyd), i.e. Q(He°)/Q(H°). (In ad- 
dition, the number of He + ionizing photons can be taken 
into account, although in case of the spectra used here, 
this contribution is small.) Figure ^ shows this ratio for 
the spectra employed in this work. The total luminosity 
of the incident continuum is set by specifying the total 
number of hydrogen ionizing photons Q(H) emitted by 
the source into An. 



In addition, this work also employs spectra appropriate to 
a whole cluster of stars. These starburst spectra were gen- 
erat ed by using the star burst synthesis code Starburst 
99 ( Leitherer et al. 1999 ). The synthesis was performed by 
adopting instantaneous star formation with Salpeter IMF 
exponent 2.35. The population of stars had an upper mass 
limit of 150M Q and a fixed stellar mass of 8.7 x 10 4 M Q . 
The evolution of the stellar cluster was described by the 
standard mass loss tracks with metallicity Z = 0.001 (i.e. 
Z/Zq — 1/20). The stellar atmospheres needed to calcu- 
late spectra were those from Kurucz (1992) and |5chmutz 
et al. 1 1992] ). A sample of spectra with different ages is 



shown in Figure |[ One may see that for longer cluster 
ages, the main contribution to the spectrum originates 
from cooler stars and the spectrum becomes softer. Note 
that in addition to the above mentioned spectra we have 
also performed test calculations with Costar stellar spec- 

as well as with spectra 



Schaerer & de Koter (199 



tra 

given by [Pauldrach et al. (1998| ). In what follows results 
for these spectra, which show generally the same trends 
as those considered in detail, are not shown. 

One parameter characterizing the "hardness" of the 
incident continuum is the ratio of the number of helium 



a -1 



T 

a -2 







"Kurucz (1991) I 


Starburst 



3.5 4 4.5 
T„ f 1 4 Kl 



5 



2 3 4 
Age [Myr] 



Fig. 4. Ratio of the number of helium ionizing photons to 
hydrogen ionizing photons as emitted by different stellar 
continuum sources. The left panel shows the single star 



spectra by Mihalas (1972 ) and Kurucz (1991 ), whereas the 
right panel refers to c ontinua of a stellar cl uster calculated 
with Starburst99 ( Leitherer et al. 1999 ) for metallicity 
Z/Z Q = 1/20 



Constraints from Observations: In this work, it is not in- 
tended to exactly model one specific cloud, but rather, 
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to simulate the whole class of low-metallicity H n-regions 
which is employed to infer helium abundances. Whether 
or not a model is suitable to describe the structure and 
emission of H n-regions used for 4 He determination may 
be decided by verifying if the physical parameters adopted 
for the Hn-regions lead to numerically calculated spec- 
tra which resemble those of the observations. Therefore, 
one may choose "realistic" models according to whether 
their emission characteristics fall broadly within certain 
observationally acceptable ranges. As a guideline for typi- 
cal line emission observed in clouds that are employed for 
prim ordial helium abun dance determinations, the sample 
from I 
tance 
to infe 
abund 
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otov ct al. (1997|) has been used. Particular impor- 
3 placed to those emission features which are used 
: main properties of the cloud as well as helium 
jjices. Ranges for these lines are defined from the 
observational sample. They cover the emission character- 
istics cf the whole sample in these specific lines. Table |] 
shows ihese ranges that the emission line ratios of a "re- 
alistic" model have to be within. 

In figure 5, each panel shows a grid of models of dif- 
ferent netallicity, where the parameters luminosity Q(H) 



Fig. 5 

space o 
ties, as 
fulfill tl 
are thof 
served 
using tl 
For 
ilar 



e 

othi :r 
Tie 



structm c 



and filling factor e are varied, respectively. The filled dots 
indicate those models, where the relative emission line 
intensities are within the ranges given in Table [l], while 
the remaining models are shown by open circles. It is evi- 
dent that the acceptable models display a band structure. 
Models along these bands are self-similar, with essentially 
identical emission spectra but varying total luminosity, in 
particular, they are described by the same Us (cf. Eq. ||). 
In the orthogonal direction (from the lower left towards 
the upper right corner) the ionization parameter Us in- 
creases. The same applies for the temperature within the 
cloud. Thus, below the band of suitable models, the line 
ratios of temperature sensitive emission lines like the for- 
bidden metal lines, tend to lower values than observed; 
whereas, models above the band show emission that is 
too strong in these lines. Towards higher metallicity, more 
models with smaller Us fulfill the range criterion. The 
emission from lower ionization stages of heavy element 
ions indicates the significance of ionization fronts at the 
outer boundaries of the cloud. These regions become more 
important relative to the main body for lower Us, as will 
be explained below. 



3. Results of the Model Calculations 

3.1. Ionization Structure 

The determination of helium abundances in H n-regions 
possibly requires a correction for unseen ionization stages 
of hydrogen or helium. When the stars ionizing the nebu- 



The luminosity Q(H) - filling factor e parameter 
"acceptable" models for two different metallici- 
labeled. Filled squares indicate the models that 
e range criterion (see Table [I]) while open circles 
e where one or more lines are not within the ob- 
i anges. All models in this figure were calculated 
Mihalas (1972] ) spectrum with T off = 45 000 K. 
spectra, the structure of the diagram stays sim- 
parameter Us stays constant along this band 
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A [A] 
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min 


max 


[On] 


3727 


0.170 


4.400 


Ne [in] .... 


3869 


0.130 


0.890 


Hel 


3889 


0.010 


0.300 


Ne [m]+H7 


3968 


0.000 


0.500 


H<5 


4101 


0.220 


0.320 


H 7 


4340 


0.410 


0.620 


[Om] 


4363 


0.030 


0.195 


Hel 


4471 


0.010 


0.060 


Hen 


4686 


0.000 


0.035 


[Om] 


4959 


0.500 


2.500 


[Om] 


5007 


1.500 


9.200 


Hei 


5876 


0.007 


0.130 


[Oi] 


6300 


0.000 


0.500 


[Shi] 


6312 


0.000 


0.037 


Ha 


6563 


2.540 


3.100 


Hei 


6678 


0.007 


0.050 


[Sn] 


6717 


0.019 


0.500 


[Sn] 


6731 


0.014 


0.374 


Hei 


7065 


0.010 


0.040 



Table 1. Range criteria for the line emission in wave- 
length A to be fulfilled by the model H n-regions. Fluxes I 
are given relative to H/3. The intervals refer to the values 



that are observed by Izotov et al. (1997). Some limits are 
relaxed by 3-10% to take into account possible deviations 
of real clouds from the idealized assumptions of the model 
calculations. 



-ionizing photons might not be sufficient to main- 1980; 


Dincrstein & Shields 1986; 


Peha 1986|; 


Armour et al. 


tdgh fraction of Hc + throughout the whole hydro- 1999; 


Vicgas et al. 200C ; 


Ballantyne et al. 200C). Recently, 



gen Stromgren sphere. Vice versa, for hard spectra of the 
ionizin g source it is also possible that a helium Stromgrcn 
sphere larger than the one of hydrogen is set up (Stasinska 



it has been attempted to co rrect for this effect , i.e. the ex- 
istence of n eutral hydrogen ( Viegas et al. 2000 ; Ballantyne 



ct al. 2000 ), though these studies reach opposite conclu- 
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sions about the induced systematic uncertainty on the in- 
ferred Y p . 

Another possible systematic error due to the ionization 
structure of H n-regions arises from the fact that the ob- 
served emission lines are always weighted towards regions 
with higher electron density. The intensity in a particular 
He I recombination line relative to a reference line like H/3 
is given by 



/(He I, A) _ / n e n Hc +a Hc {X,T)dV 
/(H/3) ~ Jn e n il +a ll (lip,T)dV 



(4) 



arc 



Here n e is the electron number density, n nc + and n K 
the number densities of He + and H + , respectively, and 
cti(A,T) are the recombination coefficients of the lines A 
considered as a function of electron temperature T. Thus, 
in regions where the electron density is small, the emission 
of a line becomes weaker even if the density of the emit- 
ting ion itself stays constant. If n e is not constant within 
the whole volume, though, the intensity ratio in Eq. (|J) 
is not proportional to He + /H + even when the electron 
temperature and, thus recombination coefficients a, stay 
constant over the cloud volume. This effect becomes par- 
ticularly important for clouds where the Stromgren sphere 
of hydrogen is smaller than the one for helium. This is due 
to the small abundances of all other elements, in partic- 
ular, H + provides the main contribution to n e . But also 
in case of a He + sphere smaller than the H + sphere this 
effect needs to be considered. The coupling of the ion- 
ization equilibrium equations of helium and hydrogen for 
photons with hv > 24.6 eV, causes variations in n e result- 



ing in smaller He I recombination line emission (3asselov 
& Gok jwirth 1995) ). 

In order to correct the observationally inferred helium 
abundances for unseen ionization stages, usually an "ion- 
ization correction factor" icf is introduced. The signifi- 
cance of this effect, is then estimated by comparing cer- 
tain emission lines that provide information on the inci- 
dent continuum and, in turn, on the ionization structure. 



In contrast to the definition of icf that is used by Viegas 
et al. (200C), Armour et al. (1999) and others, icf* as it 



is introduced here, corrects for both, unseen helium ion- 
ization stages He and He 2+ , and the effect arising from 
varying electron density within the observed volume (see 
Stasihska 1990|). The defin ition of the ionization correction 
factor byfgtasinska (1980 ), that is employed in the work of 
Pagel et al. (1992| ) and |Izotov et al. (1994 , and subsequent 
papers), takes the electron density into account. However, 
in these works icf* is assumed to be unity in most cases. 
(See the discussion later.) In the present work icf* is de- 
fined as the ratio of the true helium abundance relative 
to hydrogen and the amount of ionic helium to hydrogen 
within the volume V that is emitting line radiation at all: 



icf* 



Jn^dV /J; 



dV 



J rijidV I J n H +dV 
' Jn He +dV I Jn e n He +dV 
> / n n +dV I J n e n H +dV 



J n He dV I J n e n Rc+ dV 
J riiidV J J n e n H +dV 



(5) 



Note that the inclusion of correction for He 2+ in Eq. ((HJ) 
may be omitted, as the abundance of this ion may be 
observationally inferred from A4686 line radiation. In the 
models used here the abundance of doubly ionized helium 
is essentially negligible, such that both definitions would 
coincide. The first factor in the first equation of Eq. (JsJ) 
may be identified as the ionization correction for constant 
electron density. It is important to realize that the ef- 
fect of varying electron density, taken into account by the 
second factor in this equation, implies a required correc- 
tion which, in almost all cases, is significantly smaller (i.e. 
\{icf* — l)/(icf — 1)| < 1) than the one where the effect 
is omitted. 

3.1.1. General Trends: 

In the following, general trends of the required ionization 
correction factors for clouds with varying parameter Us 
and for ionizing sources of different spectra are studied. 
Figure || shows the relation between Us and icf* for the 
whole sample of models. The models that fulfill the range 
criterion discussed above, are shown as black filled sym- 
bols. Each track in this figure corresponds to a grid of 
models ionized by the same continuum. The lower tracks 
result from relatively hard spectra (i.e. high T eS ) while for 
softer spectra icf* tends towards larger values. (Continua 
were chosen in steps of 1000 K or 0.5 Myr, respectively.) 
Furthermore, it may be seen that when the spectrum of 
the ionizing source is fixed zc/-corrections generally be- 
come more pronounced for smaller ionization parameter 
Us- It is evident from the figure that systematic errors in- 
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Fig. 6. Ionization correction factors as a function of ion- 
ization parameter at the Stromgren sphere Us- Dark sym- 
bols refer to models fulfilling the range criteria. All models 
employed in this work are shown. 



duced by the ionization structure may, in principle, grow 
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quite large. Even for models, that show emission line char- 
acteristics similar to those observed in Hil- regions em- 
ployed for 4 He determinations (black symbols) this po- 
tential error ranges from an overestimate of 4 He by up 
to ~ 7% to an underestimate by up to > 30%. It will be 
shown below, however, that existing observational data on 
low-metallicity H n-regions seems consistent with a possi- 
ble 4 He overestimate, rather than a 4 He underestimate. 



Mihalas 




log(Q(He)/Q(H» 



Fig. 7. icf* versus Q( He°)/Q(H°) and Us for clouds ion- 
ized by [Mihalas (197^ ) spectra Dark squares refer to mod- 
els fulfilling the emission line range criterion. 



Influence of the Incident Continuum: To study the influence 
of the incident continuum on the ionization corrections in 
more detail, subsamples were analyzed. Fig. fj] shows the 
derived icf* of models ionized by different spectra as a 
function of the Q(He°)/Q(H°)-ratio and the ionization pa- 
rameter Us- Though the particular figure employs Mihalas 
spectra for the ionizing radiation, the general trends of the 
figure are hardly changed when other ionizing spectra are 
employed. This is mainly due to the choice of more phys- 
ical Q(He°)/Q(H°)-ratios as characteristics of the spec- 
tra, rather than the effective stellar temperatures. It may 
be seen that for clouds ionized by relatively soft spectra 
(small values of Q(He°)/Q(H )) the icf* is larger than 
unity for all Us ■ When the hardness of the spectrum is in- 
creased, the icf* correction approaches unity until a crit- 
ical value for Q(He°)/(3(H ) is reached. For values larger 
than about <2(He°)/(2(H ) > 0.15 required ionization cor- 
rections reverse such that values of icf* smaller than unity 
obtain. This trend is observed to be more dramatic for 
smaller Us, where the icf* may deviate significantly from 
unity. From an observational point of view it would be 
desirable to infer helium abundances within clouds with 
Us > 10~ 2 and Q(He°)/Q(H°) > 0.2 where at least within 
the simple models analyzed here icf-corrections are within 
a per cent. 



Qualitatively, these trends may be understood as fol- 
lows: For soft spectra the supply of He ionizing UV- 
photons emitted from the central source is not sufficient to 
ionize He throughout the H + sphere regardless of the ion- 
ization parameter. This leads to icf* > 1. Only photons 
with hv > 24.6 eV may ionize helium, whereas all photons 
with hv > 13.6 eV may ionize hydrogen. In the "best" 
case scenario (for maximum helium ionization of a radi- 
ation bounded nebula) all hv > 24.6 eV will be absorbed 
by helium such that the ionization equilibria of hydrogen 
and helium may be treated separately. Then for icf*> 1 
the relation 



. n Ho ( Q(He°) 



(6) 



approximately holds. This expression neglects the emis- 
sion of hydrogen ionizing radiation during helium recom- 
bination, as well as assuming Q(He°) <C Q(H°). It, never- 
theless, clarifies the existence of a critical Q(He°) /Q(H°) 
where icf* tends towards corrections larger than unity, ir- 
respective of Us- This critical value depends on the helium 
content relative to hydrogen. 

There exists an additional effect which leads to icf* 
deviating from unity, even when there is a sufficient sup- 
ply of helium ionizing photons. For decreasing Us, the 
width of the Stromgren sphere Srs increases relative to the 
Stromgren radius rs- The width Srs is given by the dis- 
tance over which ionizing radiation in partially neutral gas 
has optical depth of order unity, i.e. r w a p h e«H Srs ~ 1. 
Here a p h is an appropriate photoionization cross section. 
From this one may show that Srs/rs relates to Us via 



Srs/rs f» 5 x 10~ 7 U s 



(7) 



as long as rs 3> r$ , and under the assumption of constant 
density across the ionization front. Since the relative vol- 
ume of gas in the transition region to that in the essen- 
tially completely ionized body of the cloud is ~ 3<5rg /rs 
the effects of emission from partially ionized gas become 
increasingly important for small Us- This implies signif- 
icant ionization corrections for small Us, overestimating 
helium abundances for <5(He°)/Q(H°) > 0.15 and under- 
estimating helium abundances in the opposite case. Thus 
icf effects become asymptotically unimportant only for 
spectra that have sufficient emission of helium ionizing 
photons and for Stromgren sphere ionization parameters 
which arc sufficiently large. 

3.1.2. Observational Tools 

Because neither the ionization parameter Us itself nor the 
ratios of different Q-values are directly observable, other 
quantities have to be found in order to estimate the influ- 
ence of ionization structure on the inferred abundances. 



Radiation Softness Parameter: One frequently chosen 
method compares the abundance ratios + /0 ++ and 
S + /S ++ . This provides a measure for the quality of the 
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incident spectrum with respect to its ionizing ability. This 
"Radiation Softness Parameter" 



V 



0+ S++ 
S+ 0++ : 



(8) 



originally introduced by Vilchez fc Pagcl (1988| ), is often 
used to "read off" the ionization corr ection from mode l 
calculations like those performed by Stasihska (1990| ). 



Usually, first the particular ionic abundances are inferred 
from their respective emission lines, then rj is derived from 
these results following the definition of Eq. (||) . Here the 
particular choice of emission lines employed for the ionic 
abundance determinations may differ slightly between dif- 
ferent authors (see e.g. Pagcl ct al. 1992| ). In many obser- 
vational determinations of helium abundances, the ioniza- 
tion correction for He is simply assumed to be negligible 
for H n-regions that are ionized by sufficiently hard spec- 
tra leading to small values of rj. For instance, Pagel et al 
(1992|J] adopt no correction (i.e. icf* = 1) for those neb- 
ulae where log 77 < 0.9 and exclude objects with larger 
7/ from their analysis, [zotov fc Thuan (1998| ), more or 
less, follow this procedure. Corrections < 1 are usually 
not considered. The log 77 of the objects which this group 
uses for helium determinations, are all within the interval 
of —0.2 < log 77 < 0.4. Figure || shows the values for 77 as 
a function of the ionization parameter Us and the ratio 
Q(He°) /Q(H°). From this figure it is evident that log 77 has 
not only a dependence on the quality of the spectrum, as 
often assumed, but also on the parameter Us- H 11- regions 
which would be ideally suitable for helium abundance 
determinations would have extremely low log 77 < —0.5, 
since, following the discussion from above, they would 
have large Us and <5(He°)/Q(H°) above the critical value. 
Unfortunately, such regions seem not to be easily found in 
observational surveys. Note here that clouds illuminated 
by Kurucz (1991 ) spectra typically show larger log rj for 
the same Q(He°)/Q(H° ) and U s than those illuminated 
by eit her |Mihalas (1972[) or Starburst 99 spectra. This 



sion of photons that ionize + (hv = 2.58 Ryd) in Kurucz 
(199l|)]spectra. 

Figure|| snows the calculated icf* versus the Radiation 
Softness Parameter 77. The 77 of the H 11-regions from the 
sample of |lzotov et al. (1997 ) are within the interval that is 
indicated in this figure by the two vertical lines, the figure 
shows, that icf* may become significantly smaller than 
unity leading to a potential overestimate of the helium 
abundance in these clouds of up to ~ 6% if no correction 
is applied. 



Oxyge n-line cutoff-criterion: More recently [Ballantyne 
et al. '2000 ) suggested two interesting criteria possibly 
suitable for an estimate of the significance of ionization 
corrections. The first one involves the [O in] A5007 line rel- 
ative to H/3 and, thus, depends on the metallicity. The 
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Fig. 8. Dependency of the Radiation Softness Parameter 
77 on the ionizing photon Q(He°)/Q(H°) ratio and the ion- 
ization para meter Us ■ This f igure shows a sample of mod- 
els ionized by Mihalas (1972 ) continua of different effective 
temperatures. 
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Fig. 9. Derived values of icf* versus the Radiation 
Softness Parameter 77. The vertical lines indi cate the range 



shift is |mainly caused by the generally much smaller cmis- (1997 ) 



of 77 v alues of the sample investigated in Izotov et al 



second criterion employs the ratio [O iii]A5007/[Oi]A6300 
and should be approximately independent of metallicity: 



[Om]A5007\ 
[Oi]A6300 ) c . 



= 300 



(9) 



and 



[Oiii]A5007\ 

H /5 /cutoff 

= (0.025 ± 0.004) 



x 10 6 + (1.39 ±0.306) (10) 



They state that all objects observed with emission ratios 
lower than these cutoff values should be excluded from 
consideration since they may be subject to large "reverse" 
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ionization corrections. Figures [10| to [T^ show the results of 
our calculations in terms of these criteria. It may be seen 
that the potential error due to icf < 1 at the cutoff is still 
significant and may reach values up to 5%. From these re- 
sults it is obvious, that the suggested cutoffs are still too 
optimistic and even larger values have to be employed. 
Unfortunately, there is a scarcity of observed H n-regions 
which have sufficiently large A5007/A6300 or A5007/H/3 
to reduce «c/-corrections to a small magnitude. The main 
difference between our models and the models calculated 



1.1 



by Ballantyne et al. (2000) is the geometry. The icf cor- 
rections tend towards larger values in spherical geometry 
as the outer regions of the cloud where preferably lower 
degrees of ionization are found, contribute more to the 
total emission in a particular line than in the plane par- 
all el case. Nevertheless, o ur disagreement with the paper 
of Ballantyne et al. (2000 ) on the required cutoff values for 
negligible «c/-correction is not due to geometry. Rather, 
the proposed cutoffs by Ballantyne et al. (2000 ) are not 
s ufficient to e xclude potentially large reverse icf effects 
( [Ferland 200"o| ). 
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Fi g. 10. Mctallicity depen dent cutoff criterion suggested 
by Ballantyne et al. (200C ). Shown are the derived values 
for icf* versus the emission of [Oni]A5007 relative to H/3 
for models with low metallicity. The vertical lines indicate 
the suggested cutoff and its error. 



3.2. Temperature Variations 

In a second focus of this paper, uncertainties are investi- 
gated that may occur from the commonly used technique 
to obtain a suitable electron temperature T e for helium 
and hydrogen recombination coefficients. Even in homo- 
geneous, constant density model nebulae, the tempera- 
ture shows variations (i.e. a temperature stratification), 
whereas the abundance estimates employing emission lines 
usually assume a constant temperature. Following prior 
literature we sometimes refer to this temperature stratifi- 
cation as the existence of temperature 'fluctuations'. The 
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Fig. 11. Same as Fig. [10|but for models with high metal- 
licity. 
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Fig. 12. Mctallicity ind ependent cutoff criterion by 
Ballantyne et al. (200(1 ). This shows the icf* versus 
the line ratio [O m]A5007/[O i]A6300. The cutoff value is 
marked by the dashed line. 



variations in T e arise because the thermal equilibrium de- 
pends strongly on the local abundance of different ion- 
ization stages of coolants as, for example, oxygen. In real 
nebulae temperature variations could be even more pro- 
nounced due to, for example, more complicated density 
structures, possibly inducing 'real' small-scale fluctuations 
in the temperature. Variations of T e are typically only 
taken into account for the abundance estimate of metals, 
i.e. oxygen. This is accomplished by dividing the observed 
H II- region into two zones of different degree of ionization 
adopting different temperatures inferred from appropri- 
ate ions for each zone. The main problem in dealing with 
nebular temperatures is to find the appropriate mean tem- 
perature for the physical process considered. Since the line 
emissivities for different lines and ionization stages of dif- 
ferent elements have varying temperature dependencies, 
the bulk of the emission of certain lines used for analysis 
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may originate from either within hotter or cooler regions 
than the helium line emission. This implies that, in or- 
der to obtain accurate results for each line, a different 
temperature, determined by a appropriate weighted aver- 
age, would need to be employed (Pcimbcrt 1967). This, 



Pcimbert 



however, is a far from straightforward task (cf. 
et al. (2000)). Usually, the electron temperature influ- 
encing the emissivities of helium lines is assumed to be 
constant within the He + sphere. It is determined by ob- 
serving the ratios between the flux of collisionally excited 
oxygen lines [O in] AA5007,4959 and A4363 flAller 1984| ). 
Nevertheless, it is known that such lines are exponentially 
temperature sensitive and thus mostly originate within the 
hotter regions of the nebulae. 

In order to quantify the potential uncertainty arising 
from this method of temperature determination, another 
correction factor tcf is introduced by 

aHe(A,T[ ni]) / n e n Hc +(iU 



tcf 



_ aH(H/3,T[o m ]) fn e n H+ dV 
J n e n Hc +a Hc (X,T)dV 
J n e n H+ a H (Rp,T)dV 

Jn e n He +dV 
J n e n n +dV 



(11) 




Oin] 



where 



OH^orn]) /(He I, A) 



(12) 



N «+j Tlom] «h c (t [oi „,) i(m 

Here Tjom] denotes the "average" temperature as derived 
from the oxygen lines. The tcf factor is unity for constant 
temperature within the whole emitting volume. For "nor- 
mal" clouds this factor is smaller unity since the emission 
of [O in] lines originates in regions with high amount of 
2+ that do not coincide with the region where helium is 
ionized. Thus, this emission is weighted towards the hot- 
ter interior of the cloud. In contrast the main contribution 
to helium recombination lines is emitted within the outer, 
and thus cooler, parts of the cloud. Note, however, that in 
extreme cases the tcf may become larger than unity which 
is due to it's implicit dependence on ionization structure 
(cf. Eq. [ll]). In such cases, nevertheless, the clouds are 
already characterized by extreme icf corrections. 

Corrections due to temperature variations depend on 
the particular helium line considered since emissivities of 
different lines have different temperature dependencies. 
Observers sometimes employ the He I A6678 line alone be- 
cause it is the line that is less suspect of enhancement 
due to collisional excitation than other Hei lines. Other 
approaches use a weighted mean of all strong He I recom- 
bination lines. Of the commonly used Hei recombination 
lines A4471, A5876, and A6678, the last should be most 
sensitive to potentially arising effects from temperature 
variations. Figure [ll] shows the tcf versus r\ for all models 
computed in this work. It is derived as an equally weighted 
average of the tef's of the individual A4471, A5876, and 
A6678 lines. Because the metal lines provide the main cool- 
ing contribution, varying temperature effects are not only 



a function of ionization structure but also depend strongly 
on metallicity. Thus, the tracks in Figure ^ are broader 
than those found for the icf* corrections. From these cal- 
culation one may estimate a potential systematic overes- 
timate of Y due to a non-ideal choice of temperature of 
up to 4%. 

In case of strongly varying electron density, the effects 
of ionization structure and temperature variation are not 
anymore independent because higher order terms intro- 
duce a cross correlation between both effects. Since icf- 
effects, in some cases, may be considerably larger than tcf- 
effects this cross correlation may lead to unexpectedly high 
values for the "ideal" temperature yielding tcf > 1. Such 
models are characterized by extreme ionization structures. 
They certainly do not fulfill the simplified assumptions 
of distinct regions of high and low degree of ionization. 
Interestingly, however, even those regions show emission 
properties that may be observed in real clouds. 




Fig. 13. Temperature correction factor tcf versus the 
Radiation Softness Parameter r\ for the whole sample of 
models computed. 



In Fig.[lJ temperature correction factors are shown as 
a function of Us and the quality of the incident continuum 
as given by the ratio Q(He°)/Q(H°). The shown models 
have been computed with Mihalas spectra though results 
are very similar when different spectra are employed. It 
is evident that the tcf has it's strongest dependence on 
Us and is only mildly dependent on the spectrum. Those 
models which have tef's close to unity have large Us and 
Q(He°)/Q(H°) > 0.15 above the threshold for sufficient 
helium ionization. This is very similar to the behavior 
of ic/-effects as discussed in the last section, such that 
both potential errors could be reduced if Hn-regions of 
only very small radiation softness parameter log 77 < —0.5 
would be used for helium abundance determinations. 
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Fig. 14. tcf corrections versus the Q(He°)/Q(H°) ratio 
and Us for models ionized by Mihalas continua. 



4. Discussion 

From the results of the preceding section it is clear that 
potentially large uncertainties in the determination of he- 
lium abundances in low-metallicity H n-regions may arise 
due to the "imperfect" ionization st.mct.ure of nebulae, as 



Corrections for ionization structure seem to depend 
mostly on two parameters, the spectrum of the incident 
radiation, as well as the ratio of width to radius Srs / rs 
of the Stromgren sphere. As long as the supply of helium 
ionizing photons is sufficient (i.e. Q(He°)/<5(H°) > 0.15 
corresponding to log rj < 0.3 within the assumed mod- 
els) only "reverse" ionization corrections (icf < 1) occur. 
These ionization corrections could be, at least in principle, 
reduced to negligible levels when the ionization parameter 
at the Stromgren sphere becomes large, corresponding to 
small Srs/rg. This limit occurs for log r\ < —0.5 for the 
model H n-regions employed here. Such corrections could 
also become small when matter-bounded nebulae are con- 
sidered. In the intermediate regime —0.5 < log 77 < 0.3 icf 
corrections may become appreciable with the required cor- 
rection being larger for harder spectra. Non-uniform tem- 
perature corrections, parameterized by tcf, depend mostly 
on Srs /rs with the dependence on spectrum of secondary 
importance. These corrections also have some dependence 
on metallicity, as it is known that low-metallicity clouds 
have shallower temperature gradients. Both corrections in- 
dividually become most significant for radiation softness 
parameters log rj « leading to the apparent "triangular" 
structure in Fig. [l5| 

Ionization correction factors in helium abundance de- 



terminations have been recently investigated by Armour 



well as the existence of temperature variations Moreover 1 



both df the investigated systematic uncertainties point, 
towards a potential overestimate of helium abundances, 
since observed H 11-regions seem to be described by radia- 
tion softness parameters 77 which are indicative of possible 
"reverse" ionization corrections, rather than showing ev- 
idence .for additional amounts of neutral helium Fip; 1151 



et al (19991 ); |Viegas et al. (2000[ ) and [Ballantync et al 
(2000| ). The present work broadly confirms the poten- 



tial problem of required "reverse" ionization corrections, 
though their are differences on the viability of proposed 
criteria to exclude such uncertainties as well as the mag- 
nitude of the uncertainties when variations in electron 



shows the combined required correction factor; taking ar 

count of both effects, as a function of rj. As before, models 
which fulfill the emission line range criteria are indicated 
by black symbols. 



densit y is taken into account. Moreover, Ballantync et al 
(20001) infer an underestimate of the primordial helium 
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Fig. 15. The total required correction factor (icf* and 
tcf) versus the Radiation Softness Parameter 77 for all 
models considered in this work. 



abundance due to this effect, whereas |Vicgas et al. (2000| ) 
argue for only a moderate ^ 1% overestimate. 

Stcigman et al. (1997) concluded that the existence 
of temperature "fluctuations" in Hn-regions leads to an 
underestimate of the primordial helium abundance. This 
underestimate is claimed to be the result of an underes- 
timate of metal abundances, which leads to a decreased 
slope of helium abundance against metallicity, as well as 
to an overestimate of collisional enhancement of helium 
recombination lines. Whereas the present analysis has not 
studied the first effect, their is disagreement on the sign 
of the second effect. Steigman et al. (1997j ) lower the tem- 
peratures of H 11- regions from the inferred T[ 0n i] by a typ- 
ical amount of AT ks 2000 K, in order to account for a 
systematic overestimate of temperature from collisionally 
excited oxygen lines. An overestimate of ~ 20% in tem- 
perature leads to an overestimate of ~ 5% in the helium 
abundances since the temperature dependence of the ra- 
tio of recombination emissivities a H /a He (cf. Eq. |l2j ) is 
~ T 25 for A6678, - T 23 for A5876, and ~ T ^ for 
A4471, respectively. Coincidently, a value of ~ 5% is close 
to the largest tcf corrections found in this work. This cor- 
rcction may be compared to the potential error in derived 
|collisional enhancement factors (taken from Kingdon & 
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Fcrland (1995)) due to a possible overestimate of tem- 



perature. Most clouds in the sample of Izotov & Thuan 
(1998|)]have n e < 100 cm~ 3 . Assuming n e = 100 cm~ 3 
this change is - 0.3% for r [0lIl] 10 4 K and - 1% for 
T[om] ~ 2 x 10 4 K. Only for clouds of fairly large den- 



300 cm and temperature Tj 



[Oin] 



2 x 10 4 K 



sity n e 

does the "collisional error" become comparable ~ 2.5% 
to the error induced by employing the wrong tempera- 
ture in the recombination coefficients. It is concluded that 
an overestimate of temperature leads in most cases to an 
overestimate of helium abundance. 

Though the employed model H n-regions in this pa- 
per are probably still too simplistic to account for the 
structure of real H n-regions, it would be interesting to 
establish or refute whether existing observational data is 
consistent with systematic uncertainties of sign and mag- 
nitude as found in the previo us section. In Fig. |l7j he lium 
abundances of the sample by Izotov fc Thuan (1998 ) are 
shown as a function of the radiation softness parameter 
as inferred by these authors. It is striking to observe that 
H n-regions with log 77 > do populate the area of lower 
helium Y ~ 0.24 abundance, whereas their lower r\ coun- 
terparts seem not to. The approximate transition region of 
this trend at log rj « coincides with the peak of the trian- 
gular structure observed in Fig. |l5|, where the potentially 
largest correction factors were found. Note that the helium 
abundances shown in Fig. [l^ are not corrected for stellar 
enrichment. One may wonder if the apparent correlation 
in this graph is caused by a correlation of the inferred 
rj of H 11-regions with cloud metallicity. In Fig. |l^ radia- 
tion softness parameter is shown as a function of metal- 
licity, illustrating that such a correlation does not exist. 
An explanation of the observed trend in Fig. [l7| due to 
the existence of composite clouds, with one cloud illumi- 
nated by a hard spectrum (icf ssl) and another by a soft 
spectrum (icf > 1), seems also unlik ely, since composit e 
clouds should rarely have log rj < 0.4 ( Viegas ct al. 2000|) . 



Though it is not clear how statistically significant the 
excess of low Y H 11-regions for log rj > is, these findings 
seem suggestive that the investigated systematic uncer- 
tainties may indeed exist in the observational data. One 
may verify upon inspection of Figures |^, |^, and |IJ that 
a sample of clouds illuminated by the same hard spectra, 
but with varying Us (i.e. 5rs/rs), may not explain such 
a trend. Nevertheless, the same figures (and the knowl- 
edge that lower branches in Figure ^ correspond to larger 
Q(He°)/Q(H°), i.e. harder spectra) illustrate that a sam- 
ple of clouds with typical log Us ~ —3 illuminated by 
spectra of varying hardness may easily explain the ob- 
served trend. This is illustrated in Fig. [l^ where the com- 
bined correctio n factors are sh own for a sample of models 
illuminated by Mihalas (1972 ) and star burst spectra, and 
for two particular values of Us- If this indeed would be the 
explanation, helium abundances could be overestimated 
by a typical ~ 4%, with required corrections in individual 
cases possibly as large as ~ 8%. 
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Fig. 16. Distribution of Radiation Softness Parameter r\ 
versus metallicity in terms of the oxygen to hydrogen ratio 



for the sample by Izotov & Thuan (1998) 
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Fig. 17. Inferred helium mass fraction versus the 
Radiation Softness Parameter 77 for the sample of com- 



pact blue galaxies observed by Izotov & Thuan (1998) 



5. Conclusions 

The primordial helium abundance is commonly in- 
ferred from observationally determined helium abun- 
dances within low-metallicity extragalactic Hn-regions. 
Such observational determinations may be subject to sys- 
tematic uncertainties due to required corrections for exis- 
tence of neutral hydrogen or helium (ic/-corrections) and 
due to the possibility of non-uniform cloud temperature 
(referred to as ic/-corrections) . These effects have been in- 
vestigated in this paper. A chain of correction factors has 
been constructed that leads from the directly observable 
fluxes of helium- and hydrogen- recombination lines, and 
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to display a correlation between inferred helium abun- 
dance and radiation softness parameter in concordance 
with a typical overestimate of helium by about ~ 2 — 4%. 
In case such an interpretation of the data should prevail, 
and when uncertainties leading to a possible underesti- 
mate of helium abundances due to other effects would be 
excluded, a downward revision of inferred primordial he- 
lium would be required. 

Acknowledgements. The authors wish to acknowledge T. Abel 
and S. Buries for several helpful discussions and G. Ferland for 
communication. 



0.92 



-0.3 



-0.2 



-0.1 





log ri 



0.1 



0.2 



0.3 



Fig. 18. Resulting total correction tcf x icf* as a func- 
tion of log 7/ for models with (approximately) fixed Us, as 
labeled. Shown are models illuminated by starburst and 
Mihalas spectra. All models fulfill the range criteria. 



the inferred temperature from collisionally excited oxygen 
lines, to the true He/H ratio (cf. Eq. n3) in a nebula: 



^Hc . „ * , , ( N H + 

— icf xtcfx' 



(13) 



[Om] 



A large number of spherically symmetric model Hn- 
regions with varying parameters and illuminated by ion- 
izing radiation of different spectra have been constructed 
with the help of a photoionization code. It was shown that 
such models may be generally characterized by three phys- 
ical quantities: the number of helium- to hydrogen- ioniz- 
ing photons <5(He°)/Q(H°), metallicity, and the ionization 
parameter at the Stromgren sphere Us- It was also shown 
that Ug 1 is proportional to the ratio of width to radius 
drs/rs of the Stromgren sphere. The computed models 
are compared to the emission characteristics of observed 
H n-regions. A critical value for Q(He°)/Q(H°) > 0.15 
has been identified, such that for spectra above this crit- 
ical value only reverse icf corrections apply. Such spec- 
tra result in clouds with radiation softness parameter 
log r\ < 0.3. Reverse icf corrections may lead to a po- 
tentially large (up to ~ 6%) overestimate of helium 
abundances, particularly for Hn-regions of small log Us- 
These uncertainties may not easily be excluded even when 
[O in] A5007/[O i] A6300 -ratios ( |Ballantyne et al. 2000| ) are 
considered. In addition, a possible overestimate of helium 
abundances of similar magnitude (up to ~ 4%) may be 
present due to the existence of temperature gradients, 
provided that temperature is inferred from collisionally 
excited oxygen lines. Both effects are absent only for H n- 
regions of extraordinarily small radiation softness parame- 
ter log r\ < —0.5, corresponding to clouds of small Srs/rs 
(large Us) illuminated by sufficiently hard spectra. In con- 
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